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Abstract 

The paper aims to discuss statistical properties of the multi-agent based model of competi- 
tive growth. Each of the agents is described by growth (or decay) rule of its virtual "mass" 
with the rate affected by the interaction with other agents. The interaction depends on the 
strategy vector and mutual distance between agents and both are subjected to the agent's 
individual optimization process. Steady-state simulations yield phase diagrams with the 
high and low competition phases (HCP and LCP, respectively) separated by critical point. 
Particular focus has been made on the indicators of the power-law behavior of the mass 
distributions with respect to the critical regime. In this regime the study has revealed re- 
markable anomaly in the optimization efficiency. 
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1 Introduction 



Competitive growth [1] is one of the most generic processes observable in wide 
range of spatiotemporal scales. The rate of the growth/decay of the quantity M(t) 
which represents some "virtual mass" of the object might be defined by the first- 
order continuous dynamics 

— = Rate(M) = Growth(M) - Decay(M) . (1) 
at 
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Here the Decay (M) which contributes into overall Rate(M) comprises many pos- 
sible processes such as death, destruction, slowering due to competition. Most of 
the baseline studies start from the logistic rate form [2,3]. This minimum model 
includes the Malthusian term Growth(M) = a m M, a m > and specific quadratic 
term 

M 2 

Decay(M) = — , (2) 



where the competition is quantified by the carrying capacity [4] denoted as M c . 
Despite the logistics is well known since early studies of growth processes, for 
our purposes several comments regarding its structure will be helpful. As M(t) 
approaches M c for a m > the growth slows down until Rate(M c ) = 0. The abil- 
ity to detect or even quantify the proximity of mass saturation may be invaluable 
in a broad variety of the real world applications, enabling one to avoid the satu- 
ration. The rate reduction of the type 1 — M/M c has been used for mathematical 
description of tumor growth at tissue level [5] as well as cellular scales [2,5,6]. At 
the model level, the proximity of carrying capacity may stimulate growing entities 
to optimize their instant conditions. One can say that individual optimization [7,8] 
raises as revolutionary mechanism preventing competitors from the imposed ob- 
structions [1]. 

In the present paper we analyze statistical efficiency of optimization process of 
growing entities under the competitive conditions, where agent may change the 
strategy or migration routing to the more perspective regions by optimization. This 
goal requires much more detailed comprehensive and spatially distributed model 
than the logistic equation is, nevertheless, simple principles of logistic bound- 
ing and growth can be incorporated into construction of multiple entities called 
autonomous agents. The agent-based paradigm [9,10] states that even relatively 
simple rules may lead to very complex emergent behavior. Many examples may 
be found, but for brevity we mention only few of them: organizations of insect 
colonies [11,12], social [13], human economic behavior [10] or firms as autonomous 
entities [14]. 

The abstract agent-based model producing emergent "mass" distributions is pre- 
sented in this paper. The interpretation of M(t) or virtual "mass" itself depends 
on the scale and application field for which the model has been suggested. Here, 
the scalar can represent the variety of possible quantities, such as length or size 
of an organism, but also diameter or volume of the growing bacterial colony. We 
assume that the model can be helpful also in the study of generic features at eco- 
nomic and social scales. In this frame M(t) may represent wealth [15] of single 
seller, personal income [16], money owned by citizens [17], firm size [18]. From 
the standpoint of the total mass behavior one can, in principal, distinguish between 
the models that conserve mass [17] and models that violate this type of invariance. 
As non-conservativeness may be understood as a synonym for incompleteness or 
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missing information about all of the existing mass flows, it may be also understood 
as a realistic feature of the growth and mass-exchange models. Our focus on fluc- 
tuations is in part motivated by the effort to understand the origins of power-law 
distributions [15,16,17,18,19,20] in the systems which exhibit signs of the growth 
and competition [1,21]. 

In the economic context the interest dates back to the seminal Pareto's work [22]. 
General formulation of the problem is motivated by the fact that power-law distri- 
butions also known as Zipf's law [23,19] found also in distributions of city size 
[19]. The power-law distribution can be identified in the lifetimes [24], earthquake 
distributions [25], firm demises [14] or size of spatial colonies [26] as well. Here 
presented agent-based simulations contribute to the discussion about connection of 
power-law distributions with critical regime of driving parameters. This idea, stim- 
ulated by the concept of the feedback [27] and self-organized criticality [28,25,24], 
has been revisited in [29] . 

The structure of this paper will be as follows. In Section 2 we introduce agent- 
based model of the reactive interacting agents which manifests growth within the 
mass inequality constraints. Moreover, the agents are able to optimize their strate- 
gies and positions (see subsec. 2.3). In Section 3 statistical characteristics obtained 
by the numerical simulations are discussed. These results open the question about 
the role of the local individual optimization. Finally, the concluding remarks and 
perspectives of our approach are presented. 



2 The agent-based model 



In the paper, we incorporate mechanisms of growth, individual search and compete- 
tivness into the continuous stochastic agent-based framework and analyze equilib- 
rium statistical consequences of the complex model. Below we present model in 
more detailed focus now. 



2.1 State of agent 

The system consisting of N interacting autonomous agents, each of them equipped 
with specific abilities, is considered. It represents sufficiently complex and general 
model, where the formation and growth phenomena are interrelated with space, 
strategic and mobility issues of the competitive world. At time t, the state of ith 
agent is described by the tuple (Xf \ Sf\M^), where xf^ is the position, Sf-* the 
strategy and the mass of ith agent, respectively. The spatial coordinates are 
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taken from real space 




(3) 



The state of agent is characterized by the abstract vector of strategy 



c(*) _ r c(*) qW) cW 
°i — J i,l 5 °i,2 i ■ ■ ■ J i,d l 



SJSe<0,l>. 



(4) 



considered with normalization 




da o 



(5) 



Here S • ' defines the relative importance of the particular strategies, but the ampli- 
tude of their pursuing is determined by the mass (see in subsec. 2.2). The strategic 
vector is defined as an abstract information carrier which determines the strength 
of inter-agent interaction. In what follows, we use d x = 2 and d s = 10. 

2.2 Growth mass rules 

In analogy with logistic growth, the mass of each agent is considered to evolve 
according to discrete dynamics as it follows 



Here a is the constant growth rate parameter and (3 the feedback parameter. Gen- 
erally, the second term (— describes the effect of pairwise competition. In 
further, we consider the regime with a > 1. In that case a discrete model incorpo- 
rates the growth property. 

The competitive term in Eq.(6) is based on the overlap 



(6) 



N d B 



(7) 



i^j k=i 



weighted by the pair-wise real-space distance matrix 




(8) 
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which takes into account actual positions of agents. The tensorial structure M> ' M- 
has been chosen in analogy to scalar form of decay outlined by Eq.(2). An important 
feature of agents is their ability to perform transformation of scalar nourishment 
aMi into "vector mass" MjSj diversified by the components of Sj. In this context it 
has to be reminded that growing mass increases the impact of the respective strategy 
on the dynamics of its neighbors. Quite analogously as in the logistic function, we 
define the pairwise interaction to be proportional to the product of virtual masses 
and interaction parameter J. At the large inter-agent distances interaction turns to 
the asymptotics oc M®M® ||xf^ — X^|| -7 . For the small distances, the param- 
eter e is introduced to prevent from the proximity effects. The chosen form of the 
inter-agent pair interactions is of the short-range type (7, see in Tab. 1), spherically 
symmetric and purely repulsive. 



2.3 Optimization as a individual co evolutionary mechanism in dynamic land- 
scape 



The optimization starts by the analysis of external stimuli represented by actual fij. 
In general, individual optimization is the coevolutionary mechanism which serves 
to adapt to local competitive dynamical environment [30,9] understood as formed 
by the surrounding agents. 

In our model, the optimization of = SI (xf\ Sf } ) is considered to be the cause 
of the motion in the coordinate and strategic spaces. Here we use the hill climbing 
optimization technique which is the standard component of agent-based model- 
ing [9]. The expected consequence of the optimization is weakening of the compe- 
tition pressure on i-th agent. According to Eq.(6) smaller fij leads to slower loss 
of the mass in onward iterations. The optimization of the selected agent is applied 
with probability P opt . Technically, in the simulation process the optimization is ac- 
cepted if P opt becomes larger than random number drawn uniformly from (0, 1). 
When the step is accepted, the agent has to decide among two alternatives: posi- 
tional or strategy optimization. The individual optimization process which claims 
to find better position 

xf+^^xfUH,^)) (9) 



is accepted with probability P ps . Here H denotes the hill climbing operator charac- 
terized by N H variable displacements 5®. The iterations of strategic optimization 
formally written as 

s? +1) = #(s?UhA w ) (10) 



are accepted with probability 1 — P ps . 
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At first, Let as focus on the spatial optimization in more details. Formally, H ap- 
plied to the search for optimal coordinates Xj has been decomposed into particular 
tasks solved by the sub-operators 



Xg n) = 4ub (Xg„_i), *f>) , n = 1, 2, 3, . . . , 7V H , dD 

where X^*L = Xf ^ represents initial coordinates for the optimization process. For 
each of Nu steps we perform the calculation of corresponding for actual trial 
coordinates. It consists of calculation of J i; j matrix and overlap by means of Eq.(7) 
and Eq.(8). For kth spatial coordinate, and nth application of Eq.(ll) we suppose 
the trial move 



where Qj, , n _^ is a random number drawn from a uniform distributions (0, 1). The 
agent-based system is not conservative, but its square boundaries are impenetrable 
by mass. The agents must restrict their moves within boundaries, and thus path 
corrections are needed for xfjj^ 1 . 

The optimization uses two alternative displacements [1]: 6® G {5 X1 ,5 X2 }, where 
^xi > <5 X2 - The larger step 5 X1 is drawn with probability P big . The dichotomy of 
steps has intuitive reasons supported by the preliminary simulations. The step size 
5 Xl can be efficient to find location in the remote areas, whereas the displacement 
5 X2 helps to refine the spatial positions. 

The constant (spatially uniform) a in model does not guarantee automatically uni- 
form access to external resources as the free boundary conditions are imposed. The 
character of interactions combined with given conditions brings permanent hetero- 
geneity in the access to external sources. The nonequivalent persistent nonuniformi- 
ties at corners and edges can be in particular overcome by considering sufficiently 
large systems and short-range interactions. Formally, the optimization consists of 
the sequence of particular decisions given by 



if n(^\s?)<n(xf} n _ 1) ,s?) (is) 

then Xg n) = X ;^ ial otherwise xg n) = xg n _ 1} . 

When n = iV H the output of individual optimization X- <+1 ^ = X-*^ h ^ is achieved 
[seeEq. (9)]. 

The analogous procedure which uses two different types of steps 5 SI , 5 S , 2 is assumed 
for the optimization of Sj. The only exception is the normalization of the strategy 
vector [see Eq.(5)] which has to be applied after each optimization move. 
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The above presented agent-based model shares many features with bacterial species 
optimizing their access to nutrients in different conditions [31]. In this context the 
identification of is analogous to output of bacterial chemosensory system and the 
decision from the optimization can be the switching between staying in the same 
or moving to better place. 

2.4 Mass constraints, birth/death processes 

The iterative scheme Eq.(6) fails to describe mass dynamics when passing to ex- 
tremal values of masses. This deficiency of the model must be eliminated by ad- 
ditional constraints and limitations. We have considered the mass bounded by the 
lower cutoff, M d . If, due to competition, the mass has decreased bellow M d , the 
agent is replaced by a new one in random initial state. The situation may be inter- 
preted as death or crash. As the number of agents is conserved, the death of one 
agent opens the playground for immediate birth of his descendant with parameters 
drawn in a random way, analogously as in the stage of initialization (i.e. with initial 
mass equal to 2 M d ). The operating near threshold M d yields ergodicity gain in a 
way of extremal dynamics [28], which avoids from getting stuck. 

Preliminary numerical simulations uncovered that there remains specific issue un- 
resolved related to nonstationarity of mass distributions [32] caused by the growth 
very weighty agents. More profound analysis showed that quadratic form MiMj 
involved in Eq.(8) might not stop the growth if the masses of competing neigh- 
bours are not sufficiently high. Stationarity may be reached by introducing upper 
mass cutoff M up . By the update <— M up the agent suddenly reacts to the situ- 
ation > M up . This limitation can be understood as an extra constraint which 
guarantees renewability of sources. In the economic context the introducing of M up 
may roughly represent very restrictive taxation system. 



3 Results 

During the assembly preparation, the positions Xj and as well as vectors of strate- 
gies Sj are initialized by random values and initial values of masses is taken to 
be 2M d . Simulations have been carried out for the values of parameters listed in 
Table 1. 

One of the main purposes of the paper is to understand the impact of competition 
to the system dynamics. We decided to construct /^-dependences of the mean sta- 
tistical values. We start by rather artificial regime (5 — 0. In this trivially noncom- 
petitive situation all masses attain M up and optimization is equivalent to random 
walk. In Fig.(l)(a) we depict /9-dependent effect of the competitive reduction of the 
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Table 1 

Numerical values of model parameters used in the simulations supplemented by the short 
explanation and link to the main text where the issue is introduced. 

mean mass (M), where (...) stands for the numerical averaging over the time and 
assembly of the agents. 

More profound picture of the system behavior is achieved by analyzing the statistics 
of mass fluctuations. The mass fluctuations around the mean (M) are characterized 
by the mass dispersion = (M 2 ) — (M) 2 [see Fig. 2(a)]. The dependence of 
<7m exhibits extreme at/3 = (3 C = 1.2 x 10" 5 . The value may be interpreted as the 
critical point of the phase transition driven by (3. This scenario is consistent with the 
behaviour of the corresponding derivatives shown in insets of Fig.(l). The critical 
point separates phases of different competition (HCP and LCP, respectively). The 
phases were named according behavior of [Fig. 1(b)] . For higher f3 we observed 
better (better in a sense of smaller (fi)) way to minimize f2j that agents of HCP 
become self-improved because of more carefully chosen strategies. 

The power-law distributions are assumed to be the hallmarks of criticality [33]. 
It is therefore instructive to identify the specific features of distributions regard- 
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ing position of critical point. The mass distributions have been studied for three 
representative f3 values: (3 = 0.1/3 C , f3 c , 10/5 C [see Fig. (3)]. Despite significant 
distance from the critical regime we see that distributions may be roughly charac- 
terized by the power-law distributions with anomalies concentrated around the tails 
of high and low mass regions. The parametric robustness of the power-law like dis- 
tributions may be considered as realistic feature of the model. It is reasonable to 
suppose that anomalies originate prevailingly from the mass constraints (see sub- 
sec. 2.4). In our relatively small-size system the transition is not sharp but broad 
in (3. Also performed detailed analysis of exponents ensures parametric robustness. 
To characterize the distributions more quantitatively, the respective M dependence 
of effective local index u e s(M, 8M) has been calculated for each of them. These 
dependences have been constructed by fitting of power-law functions of some ex- 
ponent and amplitude on (— 5M + M, M + 5M) for varying central M. Any plateau 
of v cS (M, 5M) indicates the existence of particular power-law behavior on the do- 
mains larger than 8M. It should be noted that our idea of usage of local fits has been 
motivated by the phenomenology [34,35]. From the results obtained for different 
(3 we see that formation of power-law ditribution exhibits strong robustness with 
respect to (3 C in the sense that particular local intervals of the power-law behavior 
stay sufficiently far from (3 C . On the other hand, the classical concept of the phase 
transition expresses itself through minimum spread of u cS ~ 0.19 — 0.31 at the 
critical point. 

In the following we supplement our optimization results with the information about 
transport properties. The diffusion coefficient is an indirect characteristics which is 
related to the spatial individual optimization. Its definition 

1 N 

^ = 4E(Hxf +r) -Xf ) |r) W (14) 
JV i=l 

captures space and time averaging of pathways. The index (t) used to highlight 
data has been averaged over the time scale on the sub-assembly of the agents which 
live longer than r. The above mentioned averages are depicted in Fig.(l)(d). They 
indicate that HCP agents spread faster owing to reduced interactions, caused by 
their reduced mass, abating blocking from the rivals. 

The short-time efficiency of the optimization processes is analyzed using the fit- 
ness [36] function defined by the difference 

aqS*> = nfx®^, s?>) - fi(xg 0) , s?>) . (15) 

(The same differences have been calculated and averaged when the optimization of 
strategy vectors is considered.) According to the above measure, the optimization 
of the position or strategy is more efficient if the difference becomes higher. But 
the difference measure tells little about the long-time perspectives of the agents. 
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Namely at HCP the landscape varies very rapidly and optimization yields only the 
short-time benefits from the movement. In the conditions of critical regime [see 
Fig. (2) (b)] the mean (Afi) shows minimum. This confirms known fact that critical 
landscapes are landscapes of the extremal complexity [37,9]. 

The life expectancy characteristics depicted in Fig.(l)(c) indicates that the mean 
lifespan of HCP agents is considerably shorter in comparison with LCP agents. 
The combination of the above facts clearly indicates how inefficient the short-time 
optimization without prediction could be, namely in HCP constrictive conditions. 
The effect of enhanced diffusion as well as the effect of shortened life expectancy 
of agents could be considered as a typical examples of emergence, which has not 
been a priory integrated into the architecture of agents. 

Finally our study has been focused on the mean individual growth of the agent as 
a function of his lifespan. The dependences depicted in Fig. (4) demonstrate satu- 
ration in qualitative agreement with the previously mentioned logistic forms. The 
analysis shows that despite the typically nonequilibrium nature of growth models, 
due to the projection of the mass into individual lifespan coordinate (and subse- 
quent averaging over the assembly of masses belonging to the respective age), the 
mean growth dependences may be identified in the steady-state. For such numeri- 
cal output, the carrying capacity can be identified a posteriori. Since agents often 
imitate mutually destructive and self-destructive actions, it seems rather surpris- 
ing that saturation without recession appears at large lifespans. The qualitative mi- 
croexplanation is that the drop of the mass of some agents is often accompanied by 
empowered growth of their respective neighbors and thus finally low-rate growth 
stem from the most of the rivalry crowds. 
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(3 P 

Fig. 1. The mean equilibrium characteristics in LCP and HCP phases plotted as a func- 
tions of the parameter [5. The partial figures depict information about the mean: (a) agent's 
mass (M); (b) mutual overlap (fi) (match of strategies weighted by distance of agents); 
(c) life expectancy (t a ) reduced by competition at high f3; (d) diffusion coefficient D show- 
ing enhanced migration effect when the competition gets more intense. Each part of figure 
is supplemented by the inset showing how the first derivative of the corresponding quantity 
changes with respect to (3. All the anomalies are localized near the expected critical point. 
The main drawbacks here are the differences caused by the finite size effects. The depen- 
dences have been obtained by averaging of 13 independent runs going from the low to high 
(3. For each run and each fixed j3 we treated record of data corresponding to 40 000 random 
visits per agent. 



(a) (b) 




le-05 2e-05 3e-05 4e-05 5e-05 le-05 2e-05 3e-05 4e-05 5e-05 



(3 (3 

Fig. 2. The study of the fluctuations and their anomalies. The parts of figure include: 
(a) mass dispersion cjm; (b) the measure of benefits (Aft) that are gained by the short-time 
optimization. The extremes serve to identify critical point [3 C . 
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Fig. 3. The cumulative probability distributions and corresponding effective exponents 
(right hand side plots). Calculated for simulation data obtained for three representative 
values of (3: (al) [3 = &/10, (a2) /? = /3 C and (a3) /3 = 10/3 C . The local properties of 
distributions characterized by the effective exponent u e s which is a function of M (M 
is always middle point of the local fit) for two different resolutions SM. For SM = 0.5 
and corresponding (3 we found that: (bl) is cS £ (0.2,0.38); (b2) u eS G (0.2,0.31); (b3) 
v c g G (0.1, 0.42). It means that the lowest spread of local effective exponent (b2) corre- 
sponds to the critical regime. 
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Fig. 4. The return to primary motivation for the model construction [see Eq.(l)]. The parts 
showing the individual mean mass growth as a function of the lifespan of agents. Each 
age group is averaged separately. Growth dependence is calculated for three representative 
values of (3 corresponding to: (a) LCP phase, ((3 = /3 C /10); (b) critical regime; (c) HCP 
phase, (/J = 10&). 
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4 Conclusion 



In the paper the agent-based model of competitive behavior with implemented pro- 
cedure of individual optimization was investigated. Our study focused on the equi- 
librium statistics and efficiency of the individual optimization. The principal finding 
is that there exists critical regime which indicates transition from the LCP to HCP 
phase and that there are significant differences in the efficiency of optimization in 
the respective phases. In HCP the low order resources are turned to the strategically 
well organized matter. The anomaly in the efficiency belongs to the complex barri- 
ers of f2j corresponding to critical point. We observed nearly power-law behavior of 
the mass distributions robust with respect to parametric choices. The focus on the 
tails of mass distributions suggests nonequivalence of central and close to boundary 
space positions, which cause nonuniform access to the external sources. 

Among the emergent phenomena, which typically accompany the agent-based sim- 
ulations, we could mention higher mobility of lighter agents and lifespan reduced 
by their motion close to lower mass region. In the future studies we plan to inves- 
tigate the impact of an extra payoffs for the optimization and mobility which may 
strengthen competetivness. 

In the paper we present results of the equilibrium simulations of the growth. As 
the equilibrium conditions are not always suitable for the growth problems, further 
perspectives of given model can be seen in nonequilibrium applications (e.g. in 
the models of metastatic growth with dissemination of malignant cells). It would 
be also interesting to combine individualized distributed parameters, e.g. those for 
decision to optimize (agent-dependent, distributed P ps>i instead of uniform P ps ), 
and analyze their impact on the mass statistics. Further perspectives can be seen 
in the application of the realistic geografic boundary conditions, space-distributed 
sources cc(x) and assortment related to the specific sources. 
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